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The magnetic field in the solar convective zone has a random small-scale component with the r.m.s. value 
substantially exceeding the strength of a regular large-scale field. For two Majorana neutrino flavors x two 
helicities in the presence of a neutrino transition magnetic moment and nonzero neutrino mixing we analize the 
displacement of the allowed (Am 2 - sin 2 2#)-parameter region reconciled for all Underground experiments with 
solar neutrinos in dependence on the r.m.s. magnetic field value b. In contrary with the RSFP scenario with a 
regular large-scale magnetic field, we find an effective production of electron antineutrinos in the Sun even for 
small neutrino mixing through the cascade conversions like z/ e £ — > v^r — * v e R. It was found that usual SMA 
and LMA MSW parameter regions maybe forbidden while opening LOW MSW as the allowed one from the 
non-observation of v £ r in the SK experiment if random magnetic fields have strengths b > 100 kG and correlation 
lengths shorter than Lo < 1000 km. 



1. INTRODUCTION 

Recent results of the SuperKamiokande (SK) 
experiment confirmed the solar neutrino deficit at 
the level less than DATA/SSM < 0.47 @. There 
are well-known MSW, VO, RSFP solutions to the 
Solar Neutrino Problem (SNP) that have differ- 
ent signatures which could be observed through 
day/night (D/N), seasonal and 11 year periodic- 
ities of the solar neutrino flux correspondingly. 

We study here new Aperiodic Spin-Flavor Con- 
version (ASFC) scenario 

1 that, on the one hand, 
is similar with the RSFP scenario because same 
neutrino transition magnetic moment is assumed 
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while, on the other hand, differs significantly 
from the RSFP case since the presence of ran- 
dom magnetic fields leads to a large amount of 
electron antineutrinos produced efficiently by a 
non-resonant way within the convective zone of 
the Sun even for small mixing in the cascade con- 
versions V eL — > VfiR — > v eR . 

For the large mixing, sin 2 29 ~ 1, antineutrinos 
are produced through last step of such conver- 
sions for both scenarios on the way to the Earth 
in the solar wind as well as through usual vacuum 
neutrino oscillations, v e l «-> v^Li without hclicity 
change. 

Non-observation of antineutrinos in the SK ex- 
periment ($p e /$„ e < 0.035 for E e > 8.3 MeV, 
H) puts the bound on neutrino mixing angles in 
the RSFP like sin 2 29 < 0.2. We checked this 
statement from our numerical code using differ- 
ent regular magnetic field profiles || and this 
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bound almost coincides with the limit sin 2 29 < 
0.25 obtained earlier in Moreover, we con- 
firm another result in j| that the moderate mix- 
ing angle region 0.1 < sin 2 29 < 0.2 for which 
u e R maybe detectable is still acceptable for the 
RSFP solution to the SNP for the lowest al- 
lowed strengths of the magnetic field, B max < 
30 kG and for the tipical mass parameter Am 2 ~ 
10- 8 eV 2 [|. 

We make here accent on much smaller mixing 
angles, sin 2 29 <C 0.1, when the RSFP mechanism 
fails to produce v e R while our ASFC scenario re- 
mains very efficient. Moreover, the allowed Am 2 - 
region occurs quite different for the case of ran- 
dom magnetic fields. 

In all four solar neutrino experiments one mea- 
sures the integral spectrum, i.e. the number of 
neutrino events per day for the SK experiment, 

tf"=E£ ^ 0) / ' Ms)* 

a i JE th {SK) 

xa a (E) < P aa {E) > dE , (1) 

or the number of neutrino events in SNU (1 SNU 
= 10~ 36 captures per atom/per sec) in the case of 
GALLEX (SAGE) and Homestake experiments, 

j JE th (Ga,Cl) 

xa^ cl (E)<P ee (E)> dE , (2) 

where is the integral flux of neutrinos of 

kind "i" (i = pp ,Be ,pep , B) assumed to be 
constant and uniform at a given distance r = t 
from the center of the Sun, Xi (E) is the normal- 
ized differential flux, f^" 1 "*^ \i(E)dE = 1; a a 
are corresponding cross sections; the thresholds 
E th (Ga,Cl) for GALLEX (SAGE) and Homes- 
take are 0.233 MeV and 0.8 MeV correspond- 
ingly while for SK at present T t h — 5.5 MeV . 
The probabilities P aa (t) — v a{t) v a{t) : where the 
subscript a equals to a — e for v e L, a — e for v e R, 
a = jj, for v^l and a — p, for v^r correspondingly, 
satisfy the unitarity condition 

P ee (t) + P ss {t) + P„ p (t) + P m {t) - 1 , (3) 

and "< . . . >" corresponds to averaging over the 
solar interior if necessary. 



There exists, however, a problem with large 
regular solar magnetic fields, essential for the 
RSFP scenario. It is commonly accepted that 
magnetic fields measured at the surface of the Sun 
are weaker than within interior of the convective 
zone where this field is supposed to be generated. 
The mean field value over the solar disc is about 
of order 1 G and in the solar spots magnetic field 
strength reaches 1 kG. 

Because sunspots are considered to be pro- 
duced from magnetic tubes transported to the 
solar surface due to the boyancy, this figure can 
be considered as a reasonable order-of-magnitude 
observational estimate for the mean magnetic 
field strength in the region of magnetic field gen- 
eration. In the solar magnetohydrodynamics Q 
one can explain such fields in a self-consistent way 
if these fields are generated by dynamo mecha- 
nism at the bottom of the convective zone (or, 
more specific, in the overshoot layer). But its 
value seems to be too low for effective neutrino 
conversions. 

The mean magnetic field is however followed 
by a small scale, random magnetic field. This 
random magnetic field is not directly traced by 
sunspots or other tracers of solar activity. This 
field propagates through convective zone and pho- 
tosphere drastically decreasing in the strength 
value with an increase of the scale. According 
to the available understanding of solar dynamo, 
the strength of the random magnetic field inside 
the convective zone is larger than the mean field 
strength. A direct observational estimation of 
the ratio between this strengthes is not available, 
however the ratio of order 50 - 100 does not seem 
impossible. At least, the ratio between the mean 
magnetic field strength and the fluctuation at the 
solar surface is estimated as 50 (see e.g. [Q). 

This is the main reason why we consider here 
an analogous to the RSFP scenario, an aperi- 
odic spin-flavour conversion (ASFC), based on 
the presence of random magnetic fields in the so- 
lar convective zone. It turns out that the ASFC is 
an additional probable way to describe the solar 
neutrino deficit in different energy regions, espe- 
cially if current and future experiments will de- 
tect electron antineutrinos from the Sun leading 
to conclusion that neutrinos are Majorana parti- 



cles. The termin "aperiodic" simply reflects the 
exponential behaviour of conversion probabilities 
in noisy media (cf. j| , ||). 

As well as for the RSFP mechanism all argu- 
ments for and against the ASFC mechanism with 
random magnetic fields remain the same ones 
that have been summarized and commented by 
Akhmedov (sec |l(J and references therein). 

2. MAGNETIC FIELD MODEL and NU- 
MERICAL SIMULATION 

The random magnetic field is considered to be 
maximal somewhere at the bottom of convective 
zone and decaying to the solar surface. To take 
into account a possibility, that the solar dynamo 
action is possible also just below the bottom of 
the convective zone (see |(|), we accept, rather 
arbitrary, that it is distributed at the radial range 
O.7i?0 — l.Oi?0, i.e. it has the same thickness as 
the convective zone, while its correlation length 
is Lq = 10 4 km that is close to the mesogranule 
size. 

We suppose also that all the volume of the con- 
vective zone is covered by a net of rectangular do- 
mains where the random magnetic field strength 
vector is constant. The magnetic field strength 
changes smoothly at the boundaries between the 
neighbour domains obeying the Maxwell equa- 
tions. Since one can not expect the strong influ- 
ence of small details in the random magnetic field 
within and near thin boundary layers between the 
domains this oversimplified model looks applica- 
ble. 

In agreement with the SSM the neutrino source 
is supposed to be located inside 
the solar core with the radius of the order of 
R v = O.li?0 = 7 x 10 4 km. For the sake of sim- 
plicity we do not consider here the spatial distri- 
bution of neutrino emissivity from a unit of a solid 

(i) 

angle of the core image and the differences in Rl 
for different neutrino kinds " i" . Different parallel 
trajectories directed to the Earth cross different 
magnetic domains because the domain size Lg (in 
the plane which is perpendicular to neutrino tra- 
jectories) is much less than the transversal size 
= R u of the full set of parallel rays, Lq <C see 
Fig. 1. The whole number of trajectories (rays) 
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Figure 1. Geometry of neutrino trajectories in 
random magnetic fields. Finite radius of v e L- 
neutrino source, the solar core radius O.li?©, is 
shown. 

with statisticaly independent magnetic fields is 
about Rl/Ll ~ 50. 

At the stage of numerical simulation of the ran- 
dom magnetic field we generate the set of random 
numbers with a given r.m.s. value and for each 
realization of random magnetic fields along each 
of 50 rays solve the Cauchi problem for the (4 x 4) 
master equation ||. 

Then we calculate the dynamical probabilities 
Paa(t), and average them in transversal plane, i.e. 
obtain the mean arithmetic probabilities as func- 
tions of mixing parameters sin 2 29 , S. We argue 
that being additive functions of the area of the 
convective zone layer (or, the same, of the number 
of rays) these probabilities with increasing area 
become self-averaged. 

3. ASYMPTOTIC SOLUTION 

To demonstrate the properties of our model we 
briefly discuss how random magnetic field influ- 
ences the small-mixing MSW solution to SNP. 
It turns out that for the SSM exponential den- 
sity profile, typical borone neutrino energies E ~ 
7-t- UMeV, and Am 2 ~ 10~ 5 eV 2 , the MSW 
resonance occurs well below the bottom of the 
convective zone. Thus we can divide the neutrino 
propagation problem and consider two successive 
stages. First, after generation in the middle of 
the Sun, neutrinos propagate in the absence of 
any magnetic fields, undergo the non-adiabatic 
(non-complete) MSW conversion v e L — > v^l and 
acquire certain nonzero values P ee end , which 
can be treated as initial conditions at the bot- 
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torn of the convective zone. For small neutrino 
mixing S2 — > the (4 x 4)-master equation Q 
then splits into two pairs of independent equa- 
tions describing correspondingly the spin-flavor 
dynamics v e L — * v^r and v^l — > v e R in noisy 
magnetic fields. In addition, once the MSW reso- 
nance point is far away from the convective zone, 
one can also omit V e and in comparison with 
c 2 6. For v e L — > Vf^R conversion this results into a 
two-component Schrodinger equation 

VeL \ _ ( -5 fib+(t) \ ( v eL 
VfiR ) ~ \ l*b-(t) S J \ v^r 

with initial conditions 

K L | 2 (0) = P ee (0), |^fl| 2 = P w (0) = 0. As 
normalized probabilities P ee (t) and Pp,p,(t) (sat- 
isfying the conservation law P ee (t) + Pjiji(t) = 
P ee (0)) are the only observables, it is convenient 
to recast Eq.(^) into an equivalent integral form 

t t! 

S(t) = 5(0) - 4/i 2 J dbi J dt 2 S{t 2 ) x 



X{M*1)M*2) +by(t 1 )by{t 2 )] X 

x cos2(5(ti - t 2 ) - [b x (h)b y (t 2 ) - 
-b x (t 2 )b y (t 1 )]sm2S(t 1 -t 2 )}, (5) 

where S(t) = 2P ee (t) ~ P ee (0) = P ee (t) - P w (t) 
is the third component of the polarization vector 
S. 

Let us trace the line of further derivation in 
a spirit of our numerical method. Dividing the 
interval of integration into a set of equal inter- 
vals of correlation length Lo, and assuming that 
possible correlations between S(t 2 ) and bi,bj un- 
der the integral are small, Sit) itself varies very 
slowly within one correlation cell, and making use 
of statistical properties of random fields, (i) dif- 
ferent transversal components within one cell are 
independent random variables, and (n) magnetic 
fields in different cells do not correlate, we can 
average Eq.(||) thus obtaining a finite difference 
analogue. Returning back to continuous version 
we get 



5(0) exp 



4 2 sin 2 5L 

r L °jsw x 



< / < P(t') > dt' 
lo 



(6) 



where we retained possible slow space dependence 
of the r.m.s. magnetic field value. 

For S <C Lq 1 and constant r.m.s. <b 2 > we 
obtain the simple <5-correlation result || , ||, 
otherwise there remains an additional stabilizing 
factor sin 2 x/x 2 , demonstrating to what extent 
the vacuum/medium neutrino oscillations within 
one correlation cell can suppress the spin-flavor 
dynamics due to the magnetic field only. 

Another important issue is the problem of tem- 
poral dependence of higher statistical moments of 
P aa . As P aa enter the Eq.(|l]) for the number of 
events one should be certain that the averaging 
procedure does not input large statistical errors, 
otherwise there will be no room for the solar neu- 
trino puzzle itself. For 5 <C Lq 1 basing on Eq. ( 
||) we evaluate the dispersion a p. Here we present 
only the final result, the details of derivation will 
be published elsewhere: 



op = v/<P 2 e (t)>-<P ee (t) > 2 
_ P ee (0) . 



< Pee(t) > 



2V2 V 
Pee(0) 



f 



,-2rt\ 



-r< 



(7) 



(8) 



where T = (4/3)^i 2 <b 2 >L Q . We see that rela- 
tive mean square deviation of P ee from its mean 
value tends with t — * 00 to its maximum asymp- 
totic value 



a P (t) 



< Pee(t) > 



i ~ 0.707, 
y/2 



t — > 00, 



(9) 



irrespectively of the initial value P ee (0). 

This estimate is true, evidently, only for one 
neutrino ray. Averaging over N independent rays 
lowers the value (^|) in \/N times. That is for 
our case of N w 50 rays we get that maximum 
relative error should not exceed approximately 
10%, thus justifying the validity of our approach. 
For smaller magnetic fields the situation is always 
better. 

To conclude this section, it is neccesary to re- 
peat that the above estimate Eq. (||) indicates 
possible danger when treating numerically the 
neutrino propagation in noisy media. Indeed, 
usually adopted one-dimensional (i.e. along one 
ray only) approximation for the (4 x 4) master 
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equation || or (2 x 2) Eq. ^ can suffer from 
large dispersion errors and one should make cer- 
tain precautions when averaging these equations 
over the random noise before numerical simula- 
tions. Otherwise, the resulting error might be 
even unpredictable. 

4. DISCUSSION and CONCLUSIONS 

In order to find the regions in the Am 2 , sin 2 29 
plane excluded from nonobservation of v e R in 
SK it is necessary to compute and plot the iso- 
lines of the ratio $>p/$>g SM . Indeed, as antineu- 
trinos were not detected in the SK experiment 
one should claim that antineutrino flux is smaller 
than the SK background <^^{E > 8.3 MeV) < 
6 x 10 4 cm~ 2 s~ 4 [[|. Deviding this inequality 
by the SSM borone neutrino flux with energies 
E > E th = 8.3 MeV, $g SM (£ > 8.3 MeV) = 
1.7 x 10 6 cm" 2 s _1 we find the bound on the av- 
eraged over cross-section and spectrum, cascade 
transition probability v e L — > v e R > 

§^{ E > 8-3 MeV) _ 
<f>% SM (E > 8.3) ~ 

= CME )Pss (E)^E )d E k q 
f™\ B (E)a p (E)dE 

or $ P /$| SM < 3.5% . Here a p (E) = 9.2 x 
10~ 42 cm 2 [(E-1.3 MeV)/10 MeV] 2 is the cross- 
section of the capture reaction v e p — > ne + . No- 
tice that the above bound is not valid for low 
energy region below the threshold of the antineu- 
trino capture by protons, E < 1.8 MeV. 

For low magnetic fields, B , b < 20 kG, both for 
regular and random ones we do not find any vio- 
lation of the bound Eq. (|l0|). However, for strong 
magnetic fields such forbidden parameter regions 
appear in different areas over 5 and sin 2 29 for 
different kind of magnetic fields (regular and ran- 
dom). In general, the influence of random mag- 
netic fields is more pronounced as compare to the 
regular ones of the same strength. 

This property of random fields was illustrated 
in our work j|] both for regular and random mag- 
netic fields. 

It follows that the more intensive the r.m.s. 
field yj (b 2 ) in the convective zone the more effec- 



tive spin-flavour conversions lead to the produc- 
tion of the right-handed P e R, ^i?-antineutrinos. 
In the case of small mixing this is also seen from 
Eq. d). 

Our results, however, show that there exists 
strong dependence not only on the r.m.s. mag- 
netic field strength but on its correlation length 
also. In particular, when for the same b rrns = 
100 kG we substitute L — 10 3 km (a granule 
size) instead of Lq = 10 4 km (a mesogranule size 
treated in numerical simulation) the small mixing 
MSW region is fully excluded via the SK bound 
Eq. @ (see Fig. 14 in @). This is easily ex- 
plained due to the maximum value of ASFC seen 
from Eq. (g) (s'mx/x — > 1). Vice versa, if the 
small-mixing MSW solution is valid and a large 
neutrino magnetic moment exists we can extract 
from neutrino data important constraints on the 
structure of random magnetic fields deep under 
surface of the Sun. 

Extrapolating the tendency above either with 
a decrease of the correlation length Lo < 10 3 km 
or with an increase of the magnetic field values 
br.m.s. > 100 kG (or for both changes) we could 
also exclude LMA MSW region or both standard 
MSW regions would be forbidden. 

However, for such case one finds LOW MSW 
with Am 2 - 10~ 7 el^ 2 , sin 2 29 ~ 0.8 becomes to 
be allowed being forbidden for larger correlation 
lengths (Lq > 10 4 km) and lower r.m.s. values 
(6 < 100 kG) (see Fig. 14 in §). 

Thus, we develop a model of neutrino spin- 
flavour conversions in the random magnetic fields 
of the solar convective zone supposing in consis- 
tence with modern MHD models of solar mag- 
netic fields that random fields are naturally much 
higher than large-scale magnetic fields created 
and supported continuously from the small- 
scale random ones || . 

It follows that if neutrinos have a large tran- 
sition magnetic moment their dynamics in the 
Sun is governed by random magnetic fields that 
, first, lead to aperiodic and rather non-resonant 
neutrino spin-flavor conversions, and second, in- 
evitably lead to production of electron antineutri- 
nos for low energy or large mass difference region. 

Analogous results were obtained in the work 
|Tl| l where the averaging of the Redfield equation 
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for the 4x4 density matrix over (^-correlated ran- 
dom magnetic fields allowed to get simple asymp- 
totical solutions for strong random fields and led 
to the same aperiodic form of the probabilities as 
in Eq.(|. 

However, for the realistic case of solar random 
magnetic fields with a finite correlation length 
the method of direct numerical simulation of the 
Schroedinger equation given in |2) is more appro- 
priate. Then the averaging of numerical solutions 
(not Schroedinger equation!) over an ensemble of 
random field configurations (see Fig. 1) gives cor- 
rect results for this case. 

If antineutrinos & e R would be found with the 
positive signal in the BOREXINO experiment JT^] 
or, in other words, a small-mixing MSW solution 
to SNP fails, this will be a strong argument in 
favour of magnetic field scenario with ASFC in 
the presence of a large neutrino transition mo- 
ment, ji ~ 10~ n /iB for the same small mixing 
angle. 

The search of bounds on fi at the level fi ~ 
10~ n /.tB in low energy v e-scattering, currently 
planning in laboratory experiments |13|j , will be 
crucial for the model considered here. 

We would like to emphasize the impor- 
tance of future low-energy neutrino experiments 
(BOREXINO, HELLAZ) which will be sensitive 
both to check the MSW scenario and the D e R- 
production through ASFP. As it was shown in a 
recent workjl4| a different slope of energy spec- 
trum profiles for different scenarios would be a 
crucial test in favour of the very mechanism pro- 
viding the solution to SNP. 
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